FACTORIZATION OF MATRICES OF QUATERNIONS 

^ ■ TERRY A. LORING 

o 

^^ , Abstract. We review known factorization results in quaternion 

^;~V matrices. Specifically, we derive the Jordan canonical form, polar 

decomposition, singular value decomposition, the QR factorization. 
We prove there is a Schur factorization for commuting matrices, 
and from this derive the spectral theorem. We do not consider 

p^ ■ algorithms, but do point to some of the numerical literature. 

Rather than work directly with matrices of quaternions, we work 
with complex matrices with a specific symmetry based on the dual 
operation. We discuss related results regarding complex matrices 

^^ ■ that are self-dual or symmetric, but perhaps not Hermitian. 
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^ . 1. The quaternionic condition 

>■ 
(^ I It is possible to prove many factorization results for matrices of qua- 

O ■ terions by deriving them from their familiar complex counterparts. The 

amount of additional work is surprisingly small. 

J. . There are more abstract factorization theorems that apply to real 

(3 I C*-algebras, along the lines of the delightful paper by Pedersen on 

factorization in (complex) C*-algebras. We are dealing here with more 
basic questions, involving only finite-dimensional linear algebra. These 
are never (hardly ever?) addressed in basis linear algebra texts, creating 

^ ■ the impression that linear algebra over the quaternions is more difficult 

than it really is. 

There are serious hazards within linear algebra over the quaternions. 
One can be lured to difficult questions of determinants and handedness 
of the spectrum, leaving perhaps victorious, but with the impression 
that all of linear algebra over the quaternions is going to be difficult. 

We mathematicians are well-advised, when upon such uneven ground, 
to seek guidance from physicists. Such guidance helped select the top- 
ics, and helped suggest notation. 
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One topic was included for reasons of elegance, the Jordan canonical 
form. For practical purposes, the Schur decomposition will generally 
suffice. We prove that as well. 

It is assumed the reader is familiar with such things as the spectral 
theorem and the functional calculus for normal matrices. It is not 
assumed that the reader knows about C*-algebras or physics, but these 
topics are mentioned incidentally. 

Let H denote the algebra of quaternions 

m= la + bi + c] + dk a, 6, c, d e m| . 

This is an algebra over M. The canonical embedding C "^^ M sending 1 
to 1 and z to i does not make HI into an algebra over C. The trouble is 
that the embedding is not central. An additional algebraic operation 
is the involution 

la + bi + cj + dkj = a — bi — cj — dk 

which satisfies several axioms, including the following. 

a e M, a; G HI =^ (ax)* = ax* 

x*x = =^ X = 

(xy)* = y*x* 

Turning to matrices, the algebra Mjv(HI) has the expected structure 
of a unital M-algebra, plus the involution 

of conjugate-transpose. We would like to know 

A* A = =^ A = 

and 

AB = I =^ BA = I. 



We will quickly prove these after we consider a representation of M^f ^ 
on C2^. 

We have an obvious representation of MAr(HI) on H^ and quickly 
notice that since left and right scalar multiplication of H^ disagree, we 
have both left-eigenvalues A G HI solving 

A\ = Av 
and right-eigenvalues /x G HI solving 

Av = vfi. 
A glance at the survey [13j reveals that many difficulties arise. 
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Ponder the situation for real matrices, in Miv(IR)- To the chagrin 
of undergraduates, it is most natural to consider the representation 
of M7v(IR) on C^. Given a real orthogonal matrix O we get the full 
picture of why it might not diagonalize in M7v(ffi) when we look at the 
complex eigenvalues. It is at this point that we are implicitly letting 
M7v(K) act on C^. 

We do well in the case of the quaternions to regard M7v(IHI) as rep- 
resented on C^^, or in more modern terms that de-emphasizes the role 
of vectors, via a certain embedding 

X:M^(H)^M2iv(C). 

This is a very old trick. For A^ > 1 it appears to have been noticed 
first by H. C. Lee [7]. 

Definition 1.1. Given two complex N-hj-N matrices A and B we set 

We generally study complex matrices in the image of x and only at 
the end of our calculations do we draw conclusions about MAr(EI). This 
is in keeping with applications in quantum mechanics, where Hilbert 
space is always complex and time-reversal symmetry will often be incor- 
porated in the conjugate linear operation T. We define T : C^^ — )■ C^^ 
by 



:i-i) 



r 



V 

w 



-w 

V 



See [9l §2.2] for details on when T is relevant to time reversal symmetry. 

A less generous description our our approach is we plan to study 
complex matrices with a certain symmetry that is useful in physics, 
and then sell the same results to pure mathematicians by re-branding 
them as theorems about matrices of quaternions. 

The operator T is relatively well behaved, despite being only conju- 
gate linear. It preserves orthogonality, and indeed 



(1.2) 
Also 



(n,rr/) = (e,?7). 



(1.3) e^n 

which is another one-line calculation. 

Lemma 1.2. The mapping x "is well- defined, is an 
morphism, is one-to-one and satisfies 

{x{Y)r = x{Yl. 



i-algehra homo- 
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Proof. Notice that every quaternion q can be written as a + 13 j with a 
and j3 in C. This makes the map well-defined. It is clearly one-to-one 
and M-linear. Notice j3j = jfi for complex number j3 and so Bj = jB. 
Therefore 



x{A + Bj)x{C + Dj 



A B 
-B A 

AC-BD_ 
-BC 



C D 
-D C 

AD + BC 
AD -BD^ 



AC 



({AC - BD) + {AD + BC) j) 
A + Bj] (C + Df 



and 
X[A + B] 



A B 

-B ~A 



X[A*-B^J]=X[[A + B3 



U 



We can describe the image of x i^i several useful ways. We will need 
several unary operations on complex matrices. For starters we need 
the transpose A^ and the pointwise conjugate A and the conjugate- 
transpose, or adjoint, A* = A = A"^. Finally, we need a twisted 
transpose that is useful in physics. This is a "generalized involution" 
X^ called the dual operation that is defined only for X in M2Ar(C). 

Definition 1.3. For A, B, C and D all complex A^-by-A^ matrices, 
define 



A B 
C D 

Alternatively, we set X^ = 

:i.4) z 



D^ -B^ 

_^T ^T 

-ZX^Z where 




-/ 



/ 





Notice that Z could be replaced by a matrix with similar properties 
and define a variation on the dual operation. Indeed, the choice of Z 
is not standardized. 

Lemma 1.4. For X in M2Af(C) the following are equivalent. 

(1) X is in the image of x; 

(2) X* = X^; 

(3) X = — T o X oT meaning X^ = — T {XT {^)) for every vector 
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Proof. Assume 



X 



Then 



X' 



A B 
-B A 



A^ ~B^ 



xK 



W AT 

Conversely X* = X^ translates into block form as 



A* C* 
B* D* 



D^ -B^ 

_QT ^T 



SO we have proven (1) <^==^ (2). 

We compute — T o X o 7" for X again in block form, and find 



ToXoT 


V 

w 


= - 




' A B' 
C D 






-w 

V 


( 
V 


' -Aw + Bw ' 
-Cw + Dv 




' -Cw + Dv 






Aw -B\ 








= 


D 
-B 


-C ' 
A 




V 

w 





so the matrix for the linear operator — T o X o 7" is (X^) . Therefore 
(2) ^ (3). D 



Definition 1.5. Let us call 



X* = X« 



the quaternionic condition and the matrix X we will call a quaternionic 
matrix. 

We now dispose of the implications 

A* A = =^ A = 

and 

AB = I =^ BA = I 

for matrices of quaternions. These are true implications for complex 
matrices, and in particular for quaternionic matrices, and therefore true 
for matrices of quaternions. 

We pause to note some axioms of the dual operation. It behaves a 
lot like the transpose. It is linear. 



{X + aYf = X + aY^ 
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which is true even for complex a. Here X and Y are any 2N-hj-2N 
complex matrices. The dual reverses multiplication, 

{XYf = Y^XK 

It undoes itself, 

{x^y = x 

and commutes with the adjoint 

[x^y = {x*f. 

Lemma 1.6. Every matrix G in M2Ar(C) can we expressed in a unique 
way as 

G = X + iY 

with X and Y being quaternionic matrices. 
Proof. Given G we set 

X = -G^* + -G 
2 2 



and 

Y = -G^* - -G. 
2 2 

D 



2. Kramers degeneracy and Schur factorization 

Eigenvalue doubling is a key feature of self-dual, self-adjoint matri- 
ces. In physics this is called the Kramers degeneracy theorem, or the 
theory of Kramers pairs O |12l |9]. This generalizes in two ways, to give 
eigenvalue doubling given the symmetry X = X'^ and conjugate-pairing 
of eigenvalues given the symmetry X* = XK 

Such a collection of paired eigenvectors will, in good situations, form 
a unitary matrix. If U is unitary matrix that satisfies the quater- 
nionic condition then it satisfies another symmetry making it symplec- 
tic, specifically U'^ZU = Z. 

Lemma 2.1. Suppose U is a unitary 2N-by-2N matrix. The following 
are equivalent: 

(1) U is symplectic; 

(2) Ul = 0; 

(3) ZU = UZ; 

(4) UoT=ToU; 

(5) // V is column j of U for j < N then column N + j of U is 

r{v). 

Proof. This follows easily from Lemma 11.41 D 
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We need to know that the group of symplectic unitary acts transi- 
tively on C". 

Lemma 2.2. //v is any unit vector in C^^ then there is a symplectic 
unitary with Uei = v. 

Proof. Let vi = v. We use (II. 2p and (II. 3p to select in order vectors 
that are an orthonormal basis for C^^, but of the form 

Vi,Tvi,V2,Tv2,. . .,VAr,TvAr. 

If we reorder to 

Vl,V2,...,VAf,Tvi,Tv2,. . .,TvAf 

we have the columns of the desired symplectic unitary. D 

The most basic result in this realm is an ugly lemma that says that 
T maps left eigenvectors of X to right eigenvectors of X", with the 
same eigenvalue. The second part of this lemma is more elegant, if less 
general. It specifices how every quaternionic matrix has a conjugate 
symmetry it its spectral decomposition. 

Lemma 2.3. Suppose X is in M2Af(C). 

(1) IfXC = XC then 

(ro*^" = A(n)*. 

(2) IfX* = X« and X^ = X^ then 

x(n) = A(ro. 

Proof. (1) Starting with 

" A B 
C D 



X 



and 



e 



V 

w 



we find that X^ = X^ translates to 

A\ + Bw = Av 
Cv + Dw = Aw 

and (Tv)* X^ = A (Tv)* translates to 



T /iT 



w'B' +vM 

so these are equivalent conditions. 

(2) follow from (1) by taking adjoints. 



Av^ 



D 
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Now we are able to extend Kramers degeneracy to a variety of situa- 
tions, starting with a block diagonalization for commuting quaternionic 
matrices. Part (2) of Theorem 12.41 appeared in jTT] . 

Theorem 2.4. Suppose Xi, . . . ,Xk in M2Ar(C) commute pairwise and 
X* = X] for all j . 

(1) There is a single symplectic unitary U so that, for all j, 



U*XjU 






with Tj upper-triangular and Sj strictly upper-triangular. 
(2) //, in addition, the Xj are normal then there is a single sym- 
plectic unitary U so that, for all j, 



U*XjU 



Dj _0_ 
D] 



with Dj diagonal. 
(3) Every symplectic unitary has determinant one. 

Proof. (1) A finite set of commuting matrices will have a common eigen- 



vector, so let V be a unit vector so that XjV 



AjV. We know then that 



so Tv is an eigenvalue for Xj with eigenvector Xj . There is a symplectic 
unitary Ui so that UiBi = v, and then UiBn+i = Tv by Lemma \2A] 
Let Yj = UlXjUi. Then 

YjBi = XjGi 

and 

Yj^N+i = Aj-eTv+i. 
This means the column 1 and column A^ + 1 are all but zeroed-out. 



K- 



Up to a non-symplectic change of basis we are looking at block-upper 
triangular matrices that commute, so the lower-right corners in that 
basis commute. This means that the 



^. 


* 





* 





A, 





Cj 





* 


A. 


* 





5, 





D, 



(2.i: 






all commute, and each must satisfy Z* 
proven the first claim. 



Z\. By induction, we have 
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For (2) we modify the proof just a little. Starting with the Xj normal, 



we find the Yj are also normal and so 



Y, 



^. 














Aj 





Q 








A. 








B. 





D, 



This, after an appropiate basis change, would be block-diagonal, and 
from this we can conclude that the matrix in (12 .ip is normal. The 
induction proceeds as before, with the stronger conclusion that Bj = 
Cj = and Aj = Dj is diagonal. 
(3) Applying (2) we find 

~ D 0_' 
D 

where D is a diagonal unitary. Therefore 

det(V^^) = det(D) det(D) = det{D)det{D) > 

and since a unitary has determinant on the unit circle, we are done. D 



W = U 



u* 



Corollary 2.5. Every matrix in MAr(]HI) is unitarily equivalent to a 
upper-triangular matrix in 1sAn{M.) that has complex numbers on the 
diagonal. 

There is an algorithm [jLj for the Schur decomposition of quaternionic 
matrices. 



Corollary 2.6. Every normal matrix in MAr(EI) is unitarily equivalent 
to a diagonal matrix in 1sAn{C). Every Hermitian matrix in MAr(EI) 
is unitarily equivalent to a diagonal matrix in MAr(M). 

Corollary 2.7. Every Hermitian self-dual matrix X in M2Ar(C) is of 
the form 



X = U 



Di 







D, 



U* 



for some symplectic unitary U and a diagonal real matrix Dj . 

These corollaries cause us to reconsider the concept of left and right 
eigenvalues in HI and focus on just those that are in C. For details on 
left eigenvalues and non-complex right eigenvalues, consult |1]. We put 
the complex right eigenvalues in a simple context with the following. 

C^ and X e C and A, 5 G M, 



Lemma 2.8. Suppose v,w 

Then 

A B 

-B ~A 



yN 



V 

w 



A 



V 

w 
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if and only if 

(^A + Sj) (v - 3w) = (v - jw) A. 

Therefore X E C is a right eigenvalue of X E M7v(IH[) if and only if X 
is an eigenvalue o/x(X). 

Proof. This a short, direct calculation. D 

3. Jordan canonical form 

Kramers degeneracy extends to the generalized eigenvectors used to 
find the Jordan canonical form. 

Lemma 3.1. Suppose X'^ = X* . If 

(X-A/)"v = and (X-A/)""^v^O 

then 

{X-XiyTv = and {X -XlY'Wv^O. 

Proof. For any Y we have 

||rv|| = ||Fv|| = \\zYzzv\\ 

proving 

||yv|| = \\Y^*rv\\. 
Since 

{(X - M)f = (X --M)" 
the result follows. D 



We see the general idea of a proof [14] of the Jordan decomposition 
for a quaternionic matrix. When we build a Jordan basis we need 
to respect T in two ways. Whatever basis we pick for the subspace 
corresponding to A with positive imaginary part, we apply T to get the 
basis for the subspace corresponding to A. When A is real we need to 
pick generalized eigenvectors in pairs. 

The Jordan form of a quaternionic matrix is not so elegant, as each 
Jordan blocks larger than 2-by-2 gets spread around to all four quad- 
rants of the matrix. We work directly with a Jordan basis and then 
make our final conclusion in terms of quaternions. 

Theorem 3.2. Suppose X'^ = X*. There is a Jordan basis for X 
consisting of pairs of the form v, Tv. 
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Proof. Let A^a denote the subspace of all generalized eigenvectors for A 
toghether with the zero vector. Recall that just treating X as a complex 
matrix, the procedure to select a Jordan basis involves selecting, for 
each A, a basis for A^^ with the following property: whenever b is in 
this basis, then {X — A/)b is either zero or back in this basis. 

For A in the spectrum with positive imaginary part we make such 
a choice and then apply T to get a set of vectors in A^^ that has the 
correct number of elements to be a basis of Nj. It will also be a 
linearly independent set since Z and conjugation both preserve linear 
independence. Since 

(X - A/)rb = r{X - A/)b 

this basis of Nj has the desired property. 

For A in the spectrum that is real we need to modify the procedure 
for selecting the basis of A^a- A common procedure selects a basis 

(3.1) ker (X - A)^ n (ker {X - Xf-^ n (im {X - A))^ 

and constructs for Nx the Jordan basis 

|(X- A)^br,fc I 1 < ^ < ^max, 0<j<r-l, k = l,...,mrj 
Since 

r{x-xyhr,k = {x-xyrK,k 

we get the desired structure if the subspaces (13. ip are T- invariant. This 
follows from the next two lemmas and the equality 

(im(X-A))^ = ker(X*-A). 

We can now assemble the bases of the various Nx to get a Jordan basis 
built from the Kramers pairs. D 

Lemma 3.3. If X^ = X* then ker(X) is T -invariant. 

Proof. This is a restatement of the A = case of Lemma [2.3( 2). D 

Lemma 3.4. // a subspace Tio is T -invariant then {Ho) is also T- 
invariant. 

Proof. If V is in ('Ho) then for every w G Ho we have 



(Tv, w) = - (Tv, TTw) = -(v, Tw) = 0. 

D 

Corollary 3.5. Suppose X G M„(H). There is an invertihle matrix 
S G M„(EI) and a complex matrix J in Jordan Form such that X = 

s-^js. 
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4. Norms 



12 



There are two operator norms to consider on X in MAr(E[), that 
induced by quaternionic Hilbert space and that induced by complex 
Hilbert space on x i.^)- They end up identicaL 

Theorem 4.1. Suppose X is in MAr(H). Then using the norms 



^j-^^- 



on H^ and 



w 



VjVj 



on C^^ we have 



llXvl 

sup —r — TT- 
VT^O ||V|| 



sup 



IxWwl 



w 



Proof. Utihzing also the norm on C we calculate the four relevant 
norms: 



A B 
[ -B A 




V 

w 


(...])( 


V 


-] 



V 

w 

V- jw 

2 



W 



I ||2 2 

V + W 



Uv + Bw\\ + lUw - Bv\ 



Av + Bw + J {Bv - Aw) 



Uv + 5w| 



i-f: -r i|2 
\Bv - AwW 



The result follows. 

5. Singular value decomposition 



D 



As in the complex case, a slick way to prove there is a singular value 
decomposition is to work out the polar decomposition and then use the 
spectral theorem on the positive part. 

In [13j it is stated that "a little more work is needed for the sin- 
gular case" when discussing the polar decomposition. The extra work 
involves padding out a quaternionic partial isometry to be a quater- 
nionic unitary (so symplectic unitary.) 

We remind the reader that [/ is a partial isometry when {U*UY = 
U*U, or equivalently UU*U = U or U*UU* = U* or {UU*)^ = UU*. 
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If we restrict the domain and range of U we find it is an isometry from 
(ker(t/))^ to {ker{U*))^. 

Lemma 5.1. Suppose U* = U'^and U is a partial isometry in M2Ar(C). 
There is a symplectic unitary W in MAr(C) so that W^ = U^ for all 
e±ker([/). 

Proof. If V is in ker(^) then by Lemma [231 Tv is also in ker(f/). Since 
V and Tv are orthogonal, we can show that ker([/) has even dimension 
2m and that is has a basis for the form 

Vi,...,v„,Tvi,...,Tvm. 

We are working in finite dimensions so the dimension of ker(^*) is also 
2m and we select for it a basis 

Wi,...,w„,Twi,...,Twm. 

We can define W to agree with U on (ker([/)) and to send Vj to Wj 
and Tvj to Twj and so get a unitary that commutes with T, which 
means it is symplectic. D 

Lemma 5.2. Suppose X* = X" in M2Ar(C). Then there is a unitary U 
and a positive semidefinite P with U* = U'^and P* = P^ and X = UP. 

Proof. Let / be a continuous function on the positive reals with /(O) = 
and /(A) = A~2 for every nonzero eigenvalue of X*X. Then let 

W = Xf{X*X). 

The usual calculations in functional calculus tell us VF is a partial isom- 

1 
etry and that X = WP for P = {X*X)'^. Working with monomials, 

polynomials and then taking limits, we can show 

(/(y))« = / (y«) 

for any positive operator Y and so 

W'^ = f (X«X*«) X^ = f {X*X) X* = w*. 

Also pi = P = P*. 

We just showed that the matrices in the minimal polar decomposition 
are quaternionic. We use Lemma [5. II to finish the argument. D 

As expected, the polar decomposition leads to a singular value de- 
composition. 

Theorem 5.3. Suppose X* = X" in M2Ar(C). There are symplectic 
unitary matrices U and V and a diagonal matrix D with nonnegative 
real entries and D^ = D* so that X = UDV . 



R 
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Proof. We take a quaternionic polar decomposition X = WP. Since P 
is positive, we apply Theorem 12.41 to get symplectic unitary matrices 
Q and V and diagonal matrix D so that P = QDV. The eigenvalues 
of P are nonnegative, so the same is true for the diagonal elements of 
D and we have the needed factorization X = {WQ)DV. D 

6. QR FACTORIZATION 

It is easy to use Lemma [2^2] to get over EI a QR factorization theorem. 
Notice that upper triangular matrices are sent by x to matrices that 
are block upper-triangular. 

Theorem 6.1. If X* = X" in M2Af(C) the there is a symplectic unitary 
Q and R of the form 

'A B ' 
-B ~A 

with A and B upper triangular. If X & yi^iM) then there is a untiary 
Q and upper triangular matrix R in Mjv(IH[) so that X = QR. 

Proof. (2) follows directly from (1), so we prove (1). 

We apply Lemma 12.21 to the first column of X and find X = QiRi 
where 

-Bi Ai 

where Ai and Bi have zeros in their first columns, except perhaps in 
the top position. As we did earlier, we can proceed with a proof by 
induction. D 

7. Self-dual matrices 

If we study matrices with X'^ = X then we are no longer working 
directly with quaternionic matrices, but as we discuss below, there is 
a connection. We discuss a Schur factorization and a structured polar 
decomposition for self-dual matrices. The latter is a bit tricky, so we 
warm up with a structured polar decomposition for symmetric complex 
matrices. 

For dealing with a single self-dual matrix, there is the efficient Paige 
/ Van Loan algorithm [5l [10] to implement the following theorem. 

Theorem 7.1. Given Xi, . . . ,Xk in M2Ar(C) that commute pairwise 
and are self-dual, there is a single symplectic unitary U so that for all 
J, 



Ri 



U*XjU 



T^ 



with Tj upper-triangular and the Cj skew-symmetric. 
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Proof. Let v be a nonzero unit vector so that Xj\ = AjV for all j. By 
Lemma 12.31 

x;(rv) = A-(rv). 

There is a symplectic unitary Ui so that Uiei = v and L^iGtv+i = Tv. 
Let Yj = U*XjUi. Then 

Yjei = AjBi 
and 

Since e^r+i is real, we take adjoint and discover 

This means the column 1 and row A^ + 1 are all but zeroed-out, 





A, 


* 


* 


* 


j ~ 






A, 



* 
A, 










s, 


* 


A 



Basic facts about block triangular matrices show that the 



Z. 






are a commuting family of matrices, and since Ui was chosen to be 



symplectic, 
the proof. 



the Z, will be self-dual. 



As simple induction now finishes 

D 



A promising numerical technique for the joint diagonalization of two 
commuting self-dual self-adjoint matrices H and K would be to form 
the normal self-dual matrix X = A + iB and apply Paige / Van Loan 
to reduce to block diagonal form. Then apply ordinary Schur decom- 
position. This technique was used in P, §9] to diagonalize matrices 
that were exactly self-dual and approximately unitary. This idea is 
mentioned in p], section 6.5. 

It is not hard to show that the minimal polar decomposition, the one 
that is unique and can involve a partial isometry, preserves in some way 
just about any symmetry thrown at it. This is because the functional 
calculus interacts well with the dual operation [8], as well as with the 
transpose. It is a bit harder to figure what happens for the maximal 
polar decomposition. By the minimal polar decomposition is meant the 
factorization that is unique and can involve a partial isometry. By the 
maximal polar decomposition is meant the factorization that involves 
a unitary. 
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We begin with the easier result about the polar decomposition of 
complex symmetric matrices. 

Theorem 7.2. If X in M2„(C) satisfies X'^ = X then there is a 
unitary U so that U^ = U and X = U \X\. 

Proof. Again chose / with /(O) = and /(A) = A~2 for every nonzero 
eigenvalue of X*X and let 

W = Xf{X*X). 

1 

As always, ly is a partial isometry and X = WP for P = {X*X) ^ = 
\X\. Now we discover 

W^ = f (X^X) X"^ = f {XX*) X = Xf {X*X) = w. 

To create a unitary U with X = U\X\ we must extend W to map 
ker(iy) to ker (W*). We can arrange f/^ = U a.s follows. Let vi, . . . , v.^ 
be an orthonormal basis of ker(VF). Then v7, . . . , v^ will be an or- 
thonormal basis of keT{W*) = keT{W) and we define V to be zero on 
(ker(W))"^ and Vvj = vj. Thus V* will be zero on {ker^W))^ and 
V*Vj = Vj, but the same can be said about V. That means V"^ = V 
and so U = W + V will be the required symmetric unitary. 

Notice there is no structure on |X|, but considered with |X*| we get 
the formula 

(7.1) \X*\ = \X\^. 

n 

For the self-dual situation, we shall see that a similar construction 
works so long as we respect Kramers degeneracy. 

Proposition 7.3. // a partial isometry W in M2Af(C) is self-dual, 
then the initial space of W will have even dimension and T will map 
the initial space isometrically onto the final space ofW. 

Proof. Applying Lemma [2.31 to the self-adjoint matrix W*W we find 

W*Wv = v =^ {W*WfTv = Tv 
=^ {WW*)Tv = Tv 
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That is, when v G (kerW^) we have Tv G (kerVF*) and ||M^*Tv|| = 
||v||. This is useful because 

(v, W*Tv) = (v, -W*Zv) 
= {-W^Zv,y) 
= {~ZWv,y) 
= (v, W*Z^r) 
= - (v, W*rv) 

which means v and W*T"v are orthogonal, and 



w*rw*rv = w*zw*z'v 

= W*ZW^Zv 

= -W*WK 
= -W*Wv 

= —V. 

If we start with a unit vector v in (ker W) then we end up with 
an orthogonal pair of unit vectors v and w with W*T"v = w and 
W*rw = -V. 

If q is a vector in (ker W) that is orthogonal to both v and w then 
W*l~c[ will be orthogonal to both W*Tv and W*Tw since T preserves 
orthogonality everywhere and W* preserves the orthogonality of vectors 
in (ker W*) . Thus we can create a basis of (ker W) out of pairs 

^r,,W*r^r^,...,^r^,W*r^r^. 

a 

We need some examples of self-dual partial isometries. Treating vec- 
tors as 2A^-by-l matrices, if we set 

V = (Tv) w* - (Tw) V* 

then this rank two (at most) matrix is self-dual since 

V^ = -Z (-Zvw* + Zwv*)^ Z 
= —Z (wv*Z — vw*Z) Z 
= Zwv* — Zvw* 

= -r(w)v* + r(v)w* 
= v. 
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If we start with v and w orthogonal, then V will be the rank- two partial 
isometry taking v to Tw and w to Tv. We record this as a lemma 
that avoids the ugly notation. 

Lemma 7.4. Suppose v and w are orthogonal unit vectors. Then the 
partial isometry from Cv + Cw to CTv + CTw that sends v to Tw 
and w to — Tv will be self- dual. 

Theorem 7.5. If X in M2n(C) satisfies X^ = X then there is a unitary 
U so that = U andX = U \X\. 

Proof. Once more W = Xf{X*X) works to create a partial isometry 
W with X = WP for P = \X\ and this time we have W^ = W. We 
need a partial isometry from ker(VF) to ker (W*) that is self-dual. The 
dimension of ker(VF) will be even, by Lemma [7.31 Moreover if 

Vi,Wi,...,V„,W„ 

is an orthogonal basis for ker(M^) then 

Tvi, Twi, . . . , Tv„, TWra 

will be an orthogonal basis for ker(Vt^*). The needed self-dual partial 
isometry V will send Vj to Tw^ and Wj to Tvj and the self-dual unitary 
we use will he U = W + V. D 

8. The odd particle causes even degeneracy 

We hope not to scare the mathematical reader with more discus- 
sion of Kramers degeneracy. What Kramers discovered was that for a 
certain systems involving a odd number of electrons, the Hamiltonian 
always had all eigenvalues with even multiplicity. 

A mathematical manifestation of this is that the tensor product of 
two dual operations is the transpose operation in disguise. In contrast 
to that, the tensor product of three dual operations is a large dual op- 
eration in disguise. Specifically if we implement an orthogonal change 
of basis, the dual operation becomes the transpose. 

This is essentially the same as the facts that is more familiar to 
mathematicians, that H ®r H ^ M^iR) and H (g)]R H (g)]R H ^ M4(H). 

In the following, we let Z^ be as in (11. 4p . the matrix that specified 
the dual operation. 

Lemma 8.1. Consider 

U = —={I®I -iZN® Zm) ■ 
V2 

For all X G M2n{C) and Y G M2m(C), 

[/* [x^ ® ytt) f/ = ([/* (X ® Y) uf . 



FACTORIZATION OF MATRICES OF QUATERNIONS 19 

Proof. Since Z]^ = —Zk we see U = U^. Also 

U*U = -{I ®I + iZM® Zm) {I®I -iZ^® Zm) 



2 

so L'^ is a unitary and U = U^^ . 
Since [/^ = —iZ^^ ® Zm we find 

u* [x^ (g)Y^)u = u (x» (g) y«) u 

= U {-iZN ® Zm) {X^ ® Y"^) {iZn ® Zm) U 
= UU'^ {X^ Y^) U^U 
= {U* {X ® Y) U)^ . 



U 



We have refrained from discussing C*-algebras, but liere they add 
clarity. What the lemma is showing is 

(M2Jv(C)®M2M(C),tt®tt) = (M2(M+JV)(C),T) . 

Lemma 8.2. For all X e M2n{C) and Y e M2m(C), 

X^ (g)Y^ = {X(g) Y)^ 
where the H on the right is taken with respect to I ® Zm- 
Proof. This is much simpler, as 

{X ®Yf = -{I® Zm) {X ® Y)^ (/ ® Zm) 

= -{I® Zm) {X^ ® Y^) (/ ® Zm) 

= X^® {-ZmY^Zm) . 

D 

Together, these lemmas show us that a tensor product of three 
(M2Ar^(C), tl) leads to something isomorphic to 

{M2(Ni+N2+Ns){'C),\\) ■ 
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